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Abstract 

The onset mechanism of superfluidity is examined by taking the case of the capillary flow of 
liquid helium 4. In the capillary flow, a substantial fall of the shear viscosity n has been observed 
in the normal phase (T\ < T < 3.7K). In this temperature region, under the strong influence of 
Bose statistics, the coherent many-body wave function grows to an intermediate size between a 
macroscopic and a microscopic one, which is different from thermal fluctuation. We consider such 
a capillary flow by including it in a general picture that includes the flow of rotating helium 4 as 
well. Using the Kramers-Kronig relation, we express 1/n in terms of the generalized susceptibility 
of the system, and obtain a formula for the shear viscosity in the vicinity of T\. Regarding bosons 
without the condensate as a non-perturbative state, we make a perturbation calculation of the 
susceptibility with respect to the repulsive interaction. With decreasing temperature from 3.7K, 
the growth of the coherent wave function gradually suppresses the shear viscosity, and makes the 
superfluid flow stable. Comparing formulas obtained to the experimental data, we estimate that 
the ratio of the superfluid density p s (T)/p defined in the mechanical sense reaches 10 -5 just above 
T\. 

PACS numbers: 67.10.Hk, 67.25.dg, 66.20.-d 
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I. INTRODUCTION 



Superfluidity was first discovered as a frictionless flow through a capillary or a narrow 
slit [1]. The simplest experiment for this phenomenon is as follows. A capillary is set up 
vertically in gravity, the upper end of which is connected to a reservoir, and the lower end 
of which is open to a helium bath Qflfl. After some liquid is filled into the reservoir, 
the liquid flows out of the reservoir through the capillary. In gravity, the level h(t) of the 
liquid in the reservoir varies with time t as h(t) = h(0)e~ at , where a is a constant that is 
inversely proportional to the coefficient of shear viscosity rj. Figure 1 shows the kinematic 
shear viscosity v(T) = rj(T)/p(T) of liquid helium 4 in the vicinity of T\ in 1 atm (p is a 
density) When liquid helium 4 approaches the superfluid phase (r) — > 0, hence a — ► oo), 
it extremely rapidly flows down through the capillary. This phenomenon instilled us with 
the notion that the frictionless flow is a prototype of superfluidity. 

We normally interpret the phenomena in liquid helium 4 using the two-fluid model, the 
essence of which is to separate the normal-fluid and superfluid part from the beginning, 
and to assume the abrupt emergence of the superfluid part at the A-point. All anomalous 
properties above T\ are regarded as being caused by thermal fluctuations. In cooling the 
Bose systems, we observe the increase of the specific heat, or the damping of the sound 
wave at \T/T\ — 1| < 10~ 2 . These anomalous properties above T\, in which all particles 
randomly move with no specific direction, are attributed to the thermal fluctuation that is 
characteristic to the Bose systems p]. The wave functions arising from thermal fluctuation 
are randomly oriented. 

In the capillary flow, we know a remarkable experimental result in connection with this 
point. In liquid helium 4 just above the A point, the kinematic shear viscosity v(T) does 
not abruptly drops to zero at T\ as in Fig.l. Rather, after reaching a maximum value 
at 3.7K, it gradually decreases with decreasing temperature, and finally drops to zero at 
the A-point. This behavior has been often attributed to thermal fluctuations, but this 
explanation is questionable. Superfluidity, observed in the capillary flow or the rotation in a 
bucket, requires a long-lived translational motion of particles with a stable specific direction. 
The characteristic time of the flow experiment is far longer than the relaxation time r t h of 
thermal fluctuations [7|. Even if the coherent wave function obeying Bose statistics arises 
as a thermal fluctuation, it will decay long before it manifests itself as a macroscopic flow. 
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FIG. 1: The temperature dependence of the kinematic shear viscosity v(T) of liquid helium 4 
obtained in Ref. 5|. 

Furthermore, the wave function arising from thermal fluctuation is randomly oriented, but 
the macroscopic flow must have a specific direction. In this respect, it is impossible that 
thermal fluctuations lead to superfluidity [$]. Similarly, the macroscopic flow showing a 
partial disappearance of v(T) at T\ < T < 3.7K needs the stable translational motion of 
atoms with a specific direction, and is qualitatively different from thermal fluctuations. 

London stressed that the total disappearance of shear viscosity is attributed to rotv = 
over the whole volume of a liquid [9]. The partial disappearance of shear viscosity at T\ < 
T < 3.7 K in Fig.l suggests the emergence of regions in which rotv = is locally realized. 
Compared to an ordinary liquid, liquid helium 4 above T\ has a 10~ 3 times smaller coefficient 
of shear viscosity. Although in the normal phase, it is already an anomalous liquid under the 
strong influence of Bose statistics. Hence, it seems natural to assume the existence of large 
but not yet macroscopic coherent wave functions [10]. If there exist such intermediate-sized 
wave functions above T\, they must affect the macroscopic flow. 

In conventional theories, we normally do not think of the large but not yet macroscopic 
wave function. This is because we use the infinite volume limit (V — > oo) in order to make a 
clear definition of the order parameter in the phase transition. In the V — > oo limit, "large 
but not yet macroscopic" is equivalent to "microscopic", and therefore the intermediate- 
sized wave function disappears from the beginning. For the system in which the distinction 
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between "microscopic" and "macroscopic" is clear, the V — > oo limit is a realistic one. For 
the Bose system at low temperature, however, the Bose statistical coherence develops to a 
macroscopic or an intermediate size, and the container and the coherent region of a liquid are 
comparable in size. The V — > oo limit is therefore a questionable assumption for considering 
the mechanical properties of the Bose system just above T\. (For the thermodynamical 
properties, we encounter no difficulty in the V — > oo limit.) 

The shear viscosity of liquid helium 4 has been subjected to considerable experimental and 
theoretical studies. These studies are mainly focused on the famous paradox at T < T\ that 
superfluid helium 4 behaves as a viscous or non-viscous fluid, depending on the experimental 
methods. For u(T) at T\ < T < 3.7 K, it is often said that it resembles v(T) of a gas rather 
than a liquid in its magnitude and in its temperature dependence. But this explanation is 
somewhat misleading. In liquid helium 4, many features associated to Bose statistics have 
been masked by the strongly interacting nature of the liquid. Well below the lambda point 
T\, various excitations of liquid helium 4 are strictly suppressed except for phonons and 
rotons. Hence, although the excitation in a liquid, these phonons and rotons are normally 
regarded as a weakly interacting dilute Bose gas Hj. For u(T) above T\, however, the 
dilute-gas picture has no grounds, because the macroscopic condensate, which is the basis 
for the above picture, has not yet developed. To explain v(T) above T\, we must deal with 
the dissipation mechanism of a liquid and the influence of Bose statistics on it. 

At present, for the accuracy of data, the early measurements of rj(T) using the capillary 
have been supersede d by new measurements using more accurate techniques such as the vi- 
brating wire 12] 13] 14]. Above the A point, the latter method gives not only quantitatively 
similar, but more precise data of T)(T). Applying the statistical analysis to these accumu- 
lated data, a precise temperature dependence of rj(T) and v(T) above T\ was obtained Q]. 
A theory worth comparing with these precise data is desired. 

This paper gives a model of the onset of superfluid flow in the capillary. Instead of as- 
suming the macroscopic condensate from the beginning, we begin with the repulsive Bose 
system with no off-diagonal-long-range order (ODLO), and show that, with decreasing tem- 
perature above T\, the coherent many-body wave function gradually grows. In Ref.[10], we 
showed that this growth reflects in the rotational properties such as the moment of inertia 
above T\. In contrast with the rotation, however, the capillary flow of a classical liquid is 
a phenomenon that is accompanied by thermal dissipation. The frictionless capillary flow 



4 



at T < T\ is superfluidity which occurs in the strongly dissipative system when it is at 
T > T\. Hence, to consider this phenomenon, we must devise a different formalism. To take 
the nature of the liquid into consideration, we will go back to fluid mechanics in the normal 
phase, and begin with the solution of the Stokes equation for the capillary flow (Poiseuille 
formula). The capillary flow is embedded into a general picture that is applicable both to 
dissipative and to non-dissipative flows using the generalized susceptibility. Applying the 
Kramers-Kronig relation to the generalized susceptibility, we relate the capillary flow to its 
non-dissipative counterpart, rotating helium 4 in a bucket, and derive the decrease of v(T) 
from the growth of the coherent wave function. When this result is combined with that in 



Ref . [lOj , we come to understand the onset of two types of superfluidity in the shear viscosity 
and in the moment of inertia from a common origin. 

This paper is organized as follows. Section. 2 examines the capillary flow in terms of the 
generalized susceptibility. Using this method, Sec.3A gives a formula for the shear viscosity 
in the vicinity of T\, and Sec.3B examines the role of Bose statistics and the repulsive 
interaction in the suppression of the shear viscosity. Section. 4 develops a microscopic model 



of the onset of the superfluid flow. By generalizing the formalism in Ref.[10l| to the time- 
dependent case, we derive the decrease of rj just above T\ from the growth of the coherent 
wave function, and discuss the stability of superfluid flow through a capillary. To describe 
the system just above T\, the mechanical superfluid density p s (T), which is defined without 
using V — > oo, is more useful than the conventional superfluid density p s (T). In Sec. 5, using 



v(T) of Fig.l, we estimate the ratio p t 



T 



I p, and compare this p s (T) / p with another p s (T) / p 



obtained by the rotation experiment [151] • Section. 6 discusses some related problems. 



II. CAPILLARY FLOW 

In an ordinary liquid flowing along x-direction, the shear viscosity causes the shear stress 
F xy between two adjacent layers at different velocities (see Fig. 2) 

F*y = 1 W (1) 

In the linear-response theory, rj of a stationary flow is given by the following two-time 
correlation function of a tensor J xy {t) = — ^2i{Pi,xPi, y / m ) 

1 f°° 

V = J dt < J xy (0)J xy (t) > . (2) 
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FIG. 2: A flow through a capillary. 

In principle, rj of a liquid, and therefore the effect of Bose statistics on rj is obtained by 
calculating an infinite series of the perturbation expansion of Eq.(2) with respect to the 
particle interaction. However, it seems to be a hopeless attempt, because the dissipation in 
a liquid is a complicated phenomenon that allows no simple approximation jit] ]. 

Fluid mechanics is a simple phenomenological theory of non-equilibrium behaviors for 
situations in which physical quantities vary slowly in space and time. In this paper, we will 
develop an intermediate theory between the microscopic and phenomenological one 
and derive some information on the shear viscosity from the correspondence between these 
two different levels of description. In an ordinary flow through a capillary with a circular 
cross-section (a radius d and a length L), the Stokes equation has a solution of the velocity 
distribution v z (r) under the pressure difference AP such as 

, , d 2 -r 2 AP 



where r is a radius in the cylindrical coordinates (Poiseuille flow [l£j]. Fig. 2). Without loss 
of generality, we may focus on a flow velocity on the axis of rotational symmetry (z-axis). 
We define a mass-flow density j = pv(r = 0), and rewrite Eq.(3) as 



l2 AP _ p 



-*-L< " = (4) 



We call a the conductivity of a viscous liquid flow through a capillary (P = Pe z ) [19] . 
From a microscopic view point, it is a long-standing problem to derive the absolute value 



of rj using Eq.(2) |20|]. In this paper, however, instead of Eq.(2), we will regard Poiseuille's 
formula Eq.(4) as a starting point. We regard Eq.(4) as a phenomenological linear-response 
relation, and compare it to the microscopic formula for the mass- flow density j like Eq.(28). 
By this comparison, we will extract only the change of rj near T\ from the total rj in Sec. 3, 
and derive this change from a microscopic model in Sec. 4. 

In a flow through a capillary set up vertically, the level of a liquid in the reservoir varies 
with time t as h(t) = h(0)e~ at (see Appendix. A). Hence, AP(t) in Eq.(4) decreases as 
AP(t) = AP(0)e- at as well. With the aid of 

1 f°° 1 



e -at 



,2 1 ™2 



^-fidw, (5) 



(tan0 = uj/a), AP(t) is decomposed into the frequency components AP(u) within a range 
of [0, u e ] where o> e is a frequency width. Let us generalize Eq.(4) to the case of j(uj)e luJt under 
the oscillatory pressure P(u)e luJt as follows 

= -,(u,)d»^H (6) 



The conductivity spectrum ct(uj) must satisfy the following sum rule 2l|] 



a(cu)dcu = f(d), (7) 







where f(d) is a conserved quantity. (The Stokes equation gives an expression of cr(cu) 
and f(d) of the classical liquid. See Appendix. B.) In Eq.(4), we must incorporate the 
effect of time variation of AP(t) into j and a, and we define the frequency-averaged a as 
a = oj^ 1 J " e a{uj)duj using a{u) in Eq.(6) . 

The change of the starting point from Eq.(2) to Eq.(4) permits us to relate different 
types of superfluidity: (a) the frictionless capillary flow, and (b) the nonclassical flow in 
a rotating bucket. The phenomenological equation (6) is appropriate to be generalized 
so as to include different types of superfluidity. Such a method originated in studies of 
electron superconductivity [22j. Between the electrical conduction and the Meissner effect 
in superconductivity, we find a relationship parallel to the relationship between the capillary 
flow and the rotational flow in superfluidity, 

The generalization of Eq.(6) is as follows. The conductivity spectrum a(u), which de- 
scribes the dissipative response of a liquid, corresponds to an imaginary part of the gener- 
alized susceptibility xd^) that describes both dissipative and non-dissipative responses of a 



liquid. We will generalize Eq.(6) so as to include not only a response in phase with AP(w) 
but a response out of phase. Hence, <r{uS) in Eq.(6) becomes a complex number o\ + iai as 
follows 

j(cj) = - KM + ia 2 (u>)} d 2 ^M, (8) 

where cr(u) in Eq.(6) is replaced by <Ti(u). The application of the pressure gradient can be 
fitted into the general picture as follows. As an external field, we assume a fictitious velocity 



v(t) that satisfies the equation of motion 23] 



dv(t) APg) 

and replace AP(u)/L in Eq.(8) with — ipuivjui) using Eq.(9). As a result, the real and 



imaginary part of <t{u) are interchanged 241 ]. 



j(uj) = p [— u<J2{uj) + iuai(u)] d v(u). (10) 

Equation (10) is our desired formula, the coefficient of which consists not only of an imagi- 
nary part iuai(u) for the dissipative flow, but of a real part — cjo^u;) for the non-dissipative 
one, thus forming the generalized susceptibility. According to the choice of v(u), Eq.(10) 
expresses different types of flow, but the generalized susceptibility has a general structure in- 
dependent of the flow type. This means that by relating &i(uj) to 02(0;), we can derive a\{uS) 
of a capillary flow from 02(0;) of a non-dissipative flow. Causality requires that fluid particles 
begins to flow only after the pressure is applied. We can use the following Kramers-Kronig 
relation for — ua2{oj) and uai(u) in Eq.(10) 

^JbLf^^L. (12) 

7T J UJ Z — UJ IA 

If one determines <Ti{uS) in the non-dissipative flow, one can obtain <j\{u) using Eq.(ll), 
hence cr{uS) in Eq.(6). 

As a non-dissipative flow, we will consider the flow in a rotating bucket. In a rotating 
bucket, a liquid makes the rigid-body rotation owing to its viscosity. The rotational velocity, 
which is used as v{u) in Eq.(10), is a product of the angular velocity Q and the radius r 
such as v d (r) = Ct x r. For this v d (r) } the dissipation function [3] 

1 



$(r) = 2i] [ dj - -e kk 5ij ) (13) 



is zero at every r, where 2ejj = dvi/dxj + dvj/dxi. This means that, except at the boundary 
to the wall, there is no frictional force within a liquid even in the normal phase, and the 
rigid-body rotation is therefore a non-dissipative flow. (On the other hand, for the Poiseuille 
flow Eq.(3), Eq.(13) is not zero at every r except for r = 0. Fluid particles in the capillary 
flow experience thermal dissipation not only at the boundary, but within the flow.) 

The flow in a rotating bucket is formulated using the generalized susceptibility of the sys- 
tem x(r,o;) [25I (see Appendix. C), which is decomposed to the longitudinal and transverse 
part (/i, v = x, y, z) 

XrM = qj fx L {q^) + (V - qj f) X T (q^). (14) 

In the flow in a rotating bucket, the influence of the bucket propagates along the radial 
direction from the wall to the center, which is perpendicular to the particle motion driven 
by rotation (Fig. 2 of Ref. 10]). Hence, the flow in a rotating bucket is a transverse response 
of the system described by the transverse susceptibility x T (g, u). In the right-hand side of 
Eq.(10), the real part of the susceptibility —puia2{ui)d 2 is expressed as 

— pd 2 ua2(uj) = lim x T (g, cj). (15) 



Now, we can express the conductivity of the capillary flow in terms of the susceptibility 
of the system. Using Eq.(15) in the right-hand side of Eq.(ll), we obtain ai(uo) for the 
capillary flow 

12 / a 2 f°° , lim^o x T {q,^) hR v 
pd ai{u ) = / duj . (16) 

7T J ijj 1 — UJ' Z 

Quantum mechanics states that, in the decay from an excited state with an energy level 
E to a ground state with E , the higher excitation energy E causes the shorter relaxation 
time t, h/r ~ \E — Eq\. In Eq.(16), the left-hand side includes the relaxation time r in 
a i = P/i^v) = P/(4Cr), whereas the right-hand side includes the excitation spectrum in 
X T {q,uj). (We used Maxwell's relation r\ = Gr [261 ] . See Appendix.D.) In this sense, Eq.(16) 
is a many-body theoretical expression of H/r ~ \E — Eq\. 
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FIG. 3: The change of the conductivity spectrum ct(uj) from (a) in a classical fluid to (b) in a 
superfluid. cr(u) is given by Eq.(B6) in Appendix. B. 

III. SHEAR VISCOSITY OF A SUPERFLUID 



A. Phenomenological argument 

In the superfluid phase, even when the pressure difference vanishes (AP = 0) in Eq.(4), 
one observes a stable flow, hence rj = 0. The essence of superfluidity is that the normal- 
fluid and superfluid part flows without any transfer of momentum from one to the other. 
Characteristic to the frictionless flow in Eq.(6) is that in addition to the normal fluid part 
<j n (uj), the conductivity spectrum <j(uj) has a sharp peak at u = 

j(cj) = - KH + AS(u)} d 2 ^M, (17) 



where A5(oj) is a simplified expression of the sharp peak at uj = with an area A 22]. 
Figure 3 schematically illustrates such a change of c(uj) when the system passes T\. The 
separation of the sharp peak from a n (u) guarantees the absence of the momentum transfer 
(except for a n (u) near u — due to the acoustic phonon). When cr(u) has a form of 
a n {uj) + A6(w), it affects the conductivity a = J^ e a(u)du, and r] in Eq.(4) is given by 

V(T) = \ -rj^ (18) 



o 



Equation (18) is explained by the following argument on the susceptibility. 
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(1) In the classical fluid, x L (<L u ) — X T (q, uj) is satisfied at small q, and one can replace 

X T (q, uj) in Eq.(16) by x L (<?) ^) for a small g [271]. Hence, in the classical fluid the conductivity 

<j\ n {uj) of the capillary flow is given by 

, 2 / /\ 2 f 00 lining x L (g,^) , 1Q * 
<7 ln (u ) = --] duJ U 2 _ ^2 • ( 19 ) 



In the superfluid phase, however, under the strong influence of Bose statistics, the con- 
dition of x L (<L u ) — X T (<?) ^ at g — > is violated (see Sec.3B and 4). Hence, one cannot 
replace x T (<?, ^) by x L (<?, ^) in Eq.(16), and must express Eq.(16) using Eq.(19) as follows 

a x {u ) = a ln (uj ) + — — / du> H - , 20 

pixd J uj z — uj' a 

in which \im q ^ [x L (q, uj) — x T (q,uj)} corresponds to a superfluid component. For the later 
use, we define a term proportional to q^ in x^v by x^ as a quantity representing the 
balance between the longitudinal and transverse response as follows 

x^u{q,uj) = Vx {q,u) + q»q v [ — 2 

= <W T (<?,^) + Xn V (q,u). (21) 

(2) The stability of a superfluid is measured by the value of x L (<L w ) — X T ( ( li UJ ) a ^ t ne 
finite uj. In general, the superfluid flow is not perfectly stable with respect to oscillating 
external perturbations. Above a certain frequency, x L (g, cu) — x T {o.i UJ ) vanishes, and the 
system behaves as a classical fluid. (In Sec. 4, we will show this using a concrete model.) 

(3) The shear viscosity near T\ is mainly determined by the value of lirng_>o[x (?j u ) ~ 
X T {q,uj)] in the vicinity of uj = 0. In the right-hand side of Eq.(20), the finiteness of 
X L (q, uj) — x T (Qy w ) & t uj ~ leads to a sharp peak at uj' = 0, because of 

(22) 







Hence, in the superfluid phase, 0i(o/) shows a peculiar peak at uj' = like that in Fig.3.(b). 
(Since this change obeys the sum rule Eq.(7), the area A of the sharp peak is equal to the 
area enclosed by a dotted and a solid curve in Fig. 3(b).) 

(4) From now, we approximate x L (<L u ) ~ X T { < ii ljJ ) by its value at uj = in Eq.(20), and 
o~i(uj) has a form such as 

a x (u) = a ln (u) + -^lim^^O) - X T (g, 0)]S(u). (23) 
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As the temperature decreases from T\, the area of the sharp peak increases. Using Eq.(23) 
in Eq.(18), one obtains the coefficient of shear viscosity 

= (7) ! ■ (24) 



' 4 Oi^^-J^V-^^ 

Here, we define the mechanical superfluid density p s (T) = lim^otx (<Zj 0) — X T {q,fy] 
(= \im q ^ Q [x^ u (q,0)q 2 /(q^q u )}), which does not always agree with the conventional thermo- 
dynamical superfluid density p s {T). (By "thermodynamical" , we imply the quantity that 
remains finite in the V — > 00 limit.) 

Since the total density p in Eq.(24) slightly increases with decreasing temperature from 
A.2K to T\ [28], the kinematic shear viscosity v{T) = r](T)/p(T) is more appropriate than 
f]{T) to describe the change of the system around T\. We obtain a formula of v{T) = 
r](T)/p(T) as follows 

"CO = rim ' (25) 

ira z uj e p 

where u n is the kinematic shear viscosity of a classical fluid, and satisfies 0\ n = p/(Ar] n ) = 
l/(4z/ n ) in Eq.(24). In general, v n (T) of the classical liquid has the property of increasing 



monotonically with decreasing temperature 29j. In Fig.l, v(T) above 3.7i^ seems to show 



this property, but its gradual fall below 3.7K cannot be explained by the classical liquid 
picture, thereby suggesting p s (T) 7^ in Eq.(25) at T\ < T < 3.7 K. 

In the free Bose system, using j^(q,T) in Eq.(C2), x M „(g, to) in Eq.(21) has a form of 

X ^ W) = ~~v\ . + e{ P )-e{ P + q) ' (26) 

where /(e(p)) = exp(/5[e(g) + S — /i]) — 1) _1 is the Bose distribution, S is a self energy 
of a boson (we ignore to and p dependence of S by assuming it small), and /1 is a chemical 
potential. At T < T\, /(e(p)) is a macroscopic number for p = 0, and nearly zero for p ^ 0. 
Hence, in the sum over p, only two terms corresponding to p = and p = — q remain. This 
means that x^g, 0) has a form of p a {T)q IJi q v /q 2 . Practically, a nonzero p s (T) and a small 
u e vanishes v(T) at T < T\ in Eq.(25). When bosons form no condensate, however, the 
sum over p in Eq.(26) is carried out by replacing it with an integral, and q~ 2 dependence 
of X^u(q,0) disappears, hence p s (T) — and v{T) remains finite. This means that, without 
the interaction between particles, the macroscopic Bose-Einstein condensation (BEC) is the 
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necessary condition for the decrease of v(T). To explain the gradual fall of v(T) just above 
T\, we must obtain — x T (g, u;) under the particle interaction. 

B. The effect of Bose statistics and repulsive interaction 

One can physically explain the fall of the shear viscosity in Eq.(25) in terms of Bose 
statistics. When a liquid flows through a capillary, it moves in the same direction but with 
a speed that varies in a perpendicular direction. For the classical liquid, Maxwell obtained a 



simple formula rj = Gr st (Maxwell's relation. See Appendix. D) 26|, where G is the modulus 
of rigidity, and r s t is a characteristic time of structural relaxations that occur under the sheer 
stress in the flow. (The reader must not confuse this r st with T t h of thermal fluctuations.) 
This relation is useful for the interpretation of liquid helium 4 as well. In the vicinity of 
T\ in liquid helium 4, no structural transition is observed in coordinate space. Hence, G 
may be a constant at the first approximation, and therefore the fall of the shear viscosity 
is attributed mainly to the decrease of r st . In view of h/r st ~ \E — E \, the decrease of r st 
suggests the increase of the excitation energy E, and it is natural to attribute it to Bose 
statistics. The relationship between the excitation energy and Bose statistics dates back to 
Feynman's argument on the scarcity of the low-energy excitation in liquid helium 4 jsoj], in 
which he explained how Bose statistics affects the many-body wave function in configuration 
space. To the shear viscosity, we will apply his explanation. 

Consider a flow in Fig. 2, in which white circles represent an initial distribution of fluid 
particles. The long thin arrows represent the displacement from white circles on a solid 
straight line to black circles on a one-point-dotted- line curve. (The influence of adjacent 
layers in a viscous flow propagates along a direction perpendicular to the particle motion. 
Hence, the excitation associated with the shear viscosity is a transverse one.) Let us assume 
that a liquid in Fig.2 is in the BEC phase, and the many-body wave function has permuta- 
tion symmetry everywhere in a capillary. At first sight, these displacements by long arrows 
seem to be a large-scale configuration change, but they are reproduced by a set of slight 
displacements by short thick arrows from the neighboring white circles in the initial distri- 
bution. In general, the transverse modulation, such as the displacement by shear viscosity, 
does not change the particle density in the large scale, and therefore, to any given particle 
after displacement, it is always possible to find such a neighboring particle in the initial 
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distribution. In Bose statistics, owing to permutation symmetry, one cannot distinguish be- 
tween two types of particles after displacement, one moved from the neighboring position by 
the short arrow, and the other moved from distant initial positions by the long arrow. Even 
if the displacement made by the long arrows is a large displacement in classical statistics, 
it is only a slight displacement by the short arrows in Bose statistics . (In contrast, the 
longitudinal modulation results in the large-scale inhomogeneity in the particle density, and 
therefore it is not always possible to find such neighboring initial positions.) 

Let us imagine this situation in 3 N- dimensional configuration space. The excited state 
made of slight displacements, which is characteristic of Bose statistics, lies in a small distance 
from the ground state in configuration space. On the other hand, the wave function of the 
excited state is orthogonal to that of the ground state in the integral over configurations. 
Since the amplitude of the ground-state wave function is uniform in principle, the amplitude 
of the excited-state wave function must spatially oscillate around zero. Accordingly, the 
wave function of the excited state must oscillate around zero within a small distance in 
configuration space. The kinetic energy of the system is determined by the 3N-dimensional 
gradient of the many-body wave function in configuration space, and therefore this steep rise 
and fall of the amplitude raises the excitation energy of this wave function. The relaxation 
from such an excited state is a rapid process with a small r st . This mechanism explains why 
Bose statistics leads to the small coefficient of shear viscosity n = Gr st . 

When the system is at high temperature, the coherent wave function has a microscopic 
size. If a long arrow in Fig.2 takes a particle out of the coherent wave function, one cannot 
regard the particle after displacement as an equivalent of the initial one. The mechanism 
below T\ does not work for the large displacement extending over two different wave func- 
tions. Hence, the relaxation time r st changes to an ordinary long r st , which is characteristic 
of the classical liquid 31]. 

When the system is just above T\, the size of the coherent wave function is not yet macro- 
scopic, but it develops to an intermediate size. In the repulsive system with high density, 
the long-distance displacements of particles takes much energy, and therefore particles in 
the low-energy excitation are likely to stay within the same wave function. The excitation 
energy is not so high as that of the long-distance displacements, but owing to Bose statistics, 
it is not so low as that of the classical liquid. Fast relaxations within these wave functions 
is the reason for the decrease of r\ in the macroscopic capillary flow. (In contrast, the ideal 
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Bose system is too simple for such a situation to be realized, in which the macroscopic 
condensate at T < T\ is a necessary condition for the decrease of r\ as in Eq.(26).) 



IV. A MODEL OF THE ONSET OF SUPERFLUID FLOW 

Let us develop a microscopic model of the onset of superffuid flow. In deriving the shear 
viscosity rj, Eq.(4) has the following advantage over Eq.(2). In the weak-coupling system (a 
gas, and a simple liquid like liquid helium 4), the particle interaction U normally enhances 



the relaxation of individual particles to local equili 



)rium positions under a given external 



force, thereby leading to small r st and r/ = Gr st 32]. If we try to formulate this tendency 
using Eq.(2), that is, to derive a decrease of rj from an increase of U in Eq.(2), there must be 
delicate cancellation of higher-order terms in the perturbation expansion of < J X y{fyJxy{t) >■ 
On the contrary, if we regard Eq.(4), in which i] appears in the denominator of a, as a linear- 
response relation, we apply the Kubo formula not to rj but to the reciprocal 1/rj. In the 
perturbation expansion of 1/rj with respect to U, an increase of U generally leads to an 
increase ofl/rj and thereby a decrease off] . Hence, one need not expect the cancellation. 
In this respect, when we use Eq.(4) as a starting point, the influence of U on rj is naturally 
built into the formalism from the beginning. 



A. The onset of nonclassical behavior 

For liquid helium 4, we use the following hamiltonian with the repulsive interaction U 

h = J2 <pH% + ^EE K-A+^pi ( u > °)> ( 27 ) 
p p,p' <? 

where $ p denotes an annihilation operator of a spinless boson. Beginning with the repul- 
sive Bose system with no ODLO, we make a perturbation expansion. Under the particle 
interaction Hi(t), Xixv{q.i u ) is derived from 



<G\T T j^x,T)j v (0,0)\G>-- 



< 0| T T j^(x, r)j u (0, 0)exp 


r re . - 

- / drHj(T) 
. Jo 


|0 > 


< 0\exp 


- \ drH^r) 
. Jo 


|0 > 



(28) 



where (3 = l/(fc^T),r = it, and j M (x, r) is an interaction representation of Eq.(C2). 
The current-current response tensor j^x, r)j u (0, 0) (a large bubble with /i and v in 
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p+q p+q 



(a) (b) (c) 

FIG. 4: When p = p' and p + q = p' + q' in (a), an exchange of particles in (b) between ( a 
large bubble) and a bubble excitation by Hi (a small one) yields (c). 

Fig. 4 (a)) is in the medium in which particles experience the repulsive interaction: Owing 
to exp{— J Hi(r)dr) in Eq.(28), scattering of particles frequently occurs in the medium, an 
example of which is illustrated by an inner small bubble with a dotted line U in Fig.4(a). As 
the system approaches T A , j^x, r)j v (0, 0) and Hj(t) in Eq.(28) get to obey Bose statistics 
strictly, that is, the particles in the large and small bubbles in Fig. 4 (a) form a coherent wave 
function as a whole jssj : When one of the two particles in the large bubble and in the small 
bubble have the same momentum (p = p'), and when the other particles in both bubbles 
have another same momentum (p + q = p' + q'), owing to Bose statistics, a graph made by 
exchanging these particles must be included in the expansion of Eq.(28). The cutting and 
reconnection of these lines between the large and small bubbles in Fig. 4(b) yields Fig. 4(c), 
in which two bubbles with the same momenta are linked by the repulsive interaction. With 
decreasing temperature, the coherent wave function grows, and such an exchange of parti- 
cles occurs many times. Furthermore, with deceasing temperature, particles with the zero 
momentum get to play a more dominant role than others. In Fig.4, this means that taking 
only processes including zero-momentum particles becomes a good approximation: In j^jy 
of Eq.(26), a bubble with p = corresponds to an excitation from the rest particle, and a 
bubble with p = —q corresponds to a decay into the rest one. Taking only such processes 
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and continuing these exchanges, one obtains 



F p (q,iu) 



+ 



F p (q, -iu) 



(29) 



4 V [l-UFp(q,iu) l-UFp{q,-iu)\' 
where the first and second term in the right-hand side corresponds to the p = and p = —q 
term in Eq.(26), respectively, and 

(exp(/3[£ - //]) - l)- 1 - (exp(/3[e(g) + E - //]) - l)- 1 



Fp(q,iu) 



-iu + e(g) 



(30) 



Let us focus on the low-energy excitations. For a small u, the first term in the bracket 
of the right-hand side of Eq.(29) is expanded with respect to u as 
Fp(q,iu) /(e(0)) - f(e(q)) 



1 - UFp(q, iu) e(q) - U[f(e(0)) - f(e(q))} 

iu 



(31) 



1 + 



+ 



liu 



e(q)-U[f(e(0))-f(e(q))] 2(e(q) - U[f(e(0)) - f(e(q)}) 2 
Hence, Eq.(29) is expanded as 



q^ l 2Fp(q,0) 
4 Vl-UFp{q,0) 



u 



e(q)(l-UFp(q,0)) 



+ 



(32) 



Since the macroscopic flow is of our interest, let us focus on Xnv{o.i u ) at a small q. Fp(q, 0) 
in Eq.(30) is a positive and monotonically decreasing function of q 2 , which approaches zero 
as q 2 — > oo. An expansion of Fp(q, 0) around q 2 = 0, Fp(0, 0) — bq 2 + • • • has a form such as 



Fp(q,0) 







4smh ^WT)-E] 
= a — bq 2 . 



1 - 







1 



2m 



+ ■ ■ 



(33) 



For a small q, 1 — UFp(q, 0) in Eq.(32) is approximated as (1 — Ua) + Ubq 2 , and we obtain 



2a 



4 V (1 - Ua) + £/6g 2 



1 - 



u 



+ ■ ■ 



(34) 



(2m)~ 2 [(l - Ua)q 2 + Ubq*] 2 

At T ^> Ta (1 — C/a > 0), Xnv(q, 0) is approximated as 2aq^q v j [4V" (1 — £/a)]. Hence, the 
mechanical superfluid density p s (T) = Am^o^V^^)^/^^ 0)] is zero. 

With decreasing temperature, the coherent wave function— y grows. The chemical 
potential /i, which reflects the size of this wave function [34|, gradually approaches S. As 
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/j, — E — > in Eq.(30), Fp(q ) iuj) gradually increases. Hence, the higher-order terms get to 
play a dominant role in Eq.(29). As p — E — > 0, C/a (= UFp(0, 0)) in Eq.(34) increases and 
finally reaches 1, that is, 

^ = 4s.n^(^< r >- S(tf "). (35) 

At this point, Xuu(q, 0) has a form of a/(2VUb) x q^q^/q 2 - Hence, the mechanical superfluid 
density p s (T) has a nonzero value a/(2VUb). Using Ua = 1 and 6 (Eq.(33)), we obtain 

I in 

A(T) = Vsinh W T)-S]|- (36) 
If we would use the V — > oo limit, Eq.(36) is zero unless /i — £ = 0. But when we 
avoid the V — > oo limit, it remains finite even when /i - E ^ 0, and we can estimate p s /p 
using experimental data above T\ (see Sec. 5). Here, we call T satisfying Eq.(35) the onset 
temperature T on of the nonclassical behavior. 

In the vicinity of T\, Eq.(35) is approximated as Uf3 — f3 2 [p(T) — £(£7)] 2 for a small 



p — E. This condition has two solutions p(T) = E(C/) ± y/UksT. The repulsive Bose system 
is generally assumed to undergo BEC as well as a free Bose gas. Hence, with decreasing 
temperature, fJ>(T) of the repulsive Bose system rapidly increases from a negative large value, 
and reach a positive E(C/) at a finite temperature, during which course the system necessarily 



passes a state satisfying /x(T) = E(C/) — y/UksT. Consequently, T on is always above T\, and 
Ps(T) has a nonzero value Eq.(36) just above T\. This means that v(T) in Eq.(25) always 
differs from v n just above T A as in Fig.l. 

At T\ < T < T on , Eq.(36) serves as an interpolation formula of p s (T). As T approaches 
T\, p s (T) in Eq.(36) approaches (m/V)[ex.p(f3[p(T) - E]) - l]" 1 . This means that at T = 
T\, p s (T\) agrees with the conventional thermodynamical superfluid density p s {T\), and it 
abruptly reaches a macroscopic number. 



B. Stability of superfluid flow 



To discuss the stability of superfluid flow, we need a simple form of Xiwio., w), which must 
be exact at a small u, and must have a reasonable asymptotic behavior for a large uj. At 
T = T on (Ua = 1), we rewrite Eq.(34) as 



q 2 AVb 



(37) 
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D _ £ (g) 4 1 ( , V 

w tanh2 /» ( T)-E]r 



where 



V 2 

The simplest and practical form of x^iq, oo) satisfying these conditions is given by 

(39) 

Using the definition of x M „ = {q^q u / 'q 2 )[x L (q^ w ) — X T (9> c^)]; we find a quantity representing 
the balance between the longitudinal and transverse responses as follows 

1 TTt I 

X L (^)-X T (q^) = ( 10) 

The stability of superfluid flow is measured by Eq.(40). As the Bose statistical coherence 
develops in Eq.(38) — £ — > 0), the superfluid flow becomes more stable with respect 

to higher-frequency external perturbations in Eq.(40) (B — > oo) 35]. 

The stability of superfluid flow is related to the repulsive interaction U. In the repulsive 
system, when a particle is dropped from a superfluid flow by external perturbations, such 
a drop decreases the kinetic energy of the system. However, since that particle does not 
move similarly to other particles, it raises the interaction energy between particles, thereby 
raising its total energy. Hence, such a drop from the flow is prevented, and the superfluid flow 
becomes stable. In the famous argument on the stability of the superfluid by Landau, the 
change of one-particle spectrum fromp 2 /(2m) to v s p by the repulsive interaction is crucial 

361 ] . We can see the mechanism corresponding to it in Eq.(40) as follows. Even if U = in 
Eq.(34), x L {q-, u0 )~X T {q-, u ) can be written in the form of Eq.(40), with sinh 2 (|/3[/i(T) — E] |/2) 
replacing sinh \/3[fi(T) — £]]. But in this case, B is simply e(g) 2 . Hence, with increasing uo, 
X L {q,uj) — x T {q-i UJ ) vanishes far more rapidly than the case of B in Eq.(38). This means 
that the superfluid flow in the case of U — is remarkably unstable with respect to external 
perturbations. The repulsive interaction U gives l/tanh 2 (|/3[/x(T) -E]|/2) to B in Eq.(38). 
Hence, x (?) UJ )~X T {q-i w ) does not easily vanish for a large u in Eq.(40), thereby making the 
superfluid stable. In this sense, Eqs.(38) and (40) are the many-body theoretical expression 
of the mechanism pointed out by Landau. The repulsive interaction plays a significant role 
both in the emergence and in the stabilization of the frictionless flow. 
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C. Conductivity of the capillary flow in the vicinity of T\ 



As discussed in Sec. 3, what exists behind the suppression of shear viscosity near T\ is 
the change of conductivity spectrum a{oj) near T\. In addition to the sharp peak at u = 0, 
there must be a gradual change in a(uo). Using Eq.(40) in Eq.(20) with the aid of 

1 1 -7T 



dcu 



uj 2 — uj' 2 



(41) 



we obtain 



(42) 



+ 



v- 1 

pud 2 



m 



sinh|/^(T)-£]| 



x , , 4 ( k B T 



7T" 



m 



7T 



cosh2 ^(T)-m 



B 



where q e in e(q) and B in Eq.(42) is a wave number corresponding to the characteristic 
length of a rotating bucket or a capillary (from now, we simply denote a\{uS) by <t{uj)). 
a\ n (uo) in the right-hand side of Eq.(42) represents the conductivity of a classical liquid 
(<j(oo) in Fig.3(a)). The sharp peak and the negative continuous part in the bracket of 
Eq.(42) represent the change of a(uo) from Fig.3(a) to 3(b). The total conductivity <j(uo) 
must satisfy the sum rule Eq.(7) regardless of whether it is in the normal or the superfluid 
phase, and therefore we introduced a normalization constant 4/-7T 2 in Eq.(42). (With the 
aid of J* °° duj/{l + u 2 '/ ' B) = ti\^B/2, the second term in the bracket of the right-hand side 
of Eq.(42) yields a term proportional to 1/ sinh \@[fj,(T) — E]| in J °° a(uj)duj, and 4/ir 2 is 
determined so that the first and second term cancel each other.) 

(Ti n (w) in Eq.(42) is given by the real part of Eq.(B5) in Appendix.B. Hence, the conduc- 
tivity spectrum a n {uj) of the normal-fluid flow has a form such as 

/ \ 

Jin I , N (43) 



a n (uj 



ujd 2 
4 



V 



n 2 nd 2 \e(q e ) 2 s ^ 



1 



00 

1+ B 



(Jo is the zero-th order Bessel function, and n = p/m). The temperature dependence 
comes from second term in the right-hand side of Eq.(43). As p(T) — S — > 0, the second 
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term slightly increases owing to 1/ cosh 2 \(3[p(T) — S]/2|, but owing to B — > oo, J °° a n (u)duj 
deceases. On the other hand, a s (u) of the superfluid flow satisfies 

oo l y-l 

< u)du = ^ sinh|/3[MT)-S]| - (44) 

The the sharp peak at uo = increases owing to 1/ sinh \(3[p(T) — Equation (43) and 
(44) describe the change of the conductivity spectrum o"(u;) in the vicinity of T\ in Fig.3. 

V. COMPARISON TO EXPERIMENTS 

Let us estimate p s (T)/p using Eq.(25) and Fig.l, and compare it to Eq.(36). 
(1) Comparing Fig.l and Eq.(25), we regard v n in the right-hand side of Eq.(25) as ^(T on ) 
with T on = 3.7K, and obtain p s (T)/p by 

* P S (T) mJfr 



ird 2 u e p 



v{TY l -v{T on )-\ (45) 



(a) In the experiments using a capillary, the level of a liquid in the reservoir varies as 
h{t) = h(0)e~ at . In Ref.[4], this type of experiment was performed in liquid helium 3, and 
a of liquid helium 3 at 1.105i^ was 5 x 10 s , which value must be close to that at 3.7K. 
Furthermore, liquid helium 3 and 4 have close values of a under the same experimental 
condition near 3.7 K (see Appendix. A.) Hence, we approximately use 5 x 10 _4 s _1 for liquid 
helium 4 at 3.7A". Using this a in Eq.(5), we estimate the frequency width u e of AP(u) in 
AP(0)e~ at to be the half width of the sharp peak 1/Vlu 2 + a 2 at u — 0, obtaining uj e ~ \/3a. 
A rough estimate of uj e is uj e ~ 5 x 10~ 4 rad/ s. (b) For the capillary radius d, we use a typical 
value d ~ 10~ 2 cm in Ref.[4|. (c) In Fig.l, i/(T) _1 — z/(T on ) _1 reaches 2 x 10 3 s/cm 2 just above 
T A . Using these u e and d in comparing Eq.(45) to v{T) of Fig.l, we find just above T A 

I — I 

Figure. 5 shows the temperature dependence of pjT)/p given by v{T) of Fig.l [371 ] . 



In the rotation experiment by Hess and Fairbank 15j, the moment of inertia I z just above 



T\ is slightly smaller than the normal phase value If (see Sec. 6. A). Using these currently 
available data, Ref. jlo| roughly estimates p s (T)/p as p s (T\ + 0.03K)/p = 8 x 10~ 5 , and 
p s (T\ + 0.28K)/p = 3 x 10 -5 . Although the precision of these estimates is limited, the 
agreement of p s (T)/p in Fig.5 with these values is good beyond our expectation. 
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FIG. 5: p s {T)lp obtained by Eq.(45) using u(T) in Fig.l. 

(2) In Sec. 4, we obtained the interpolation formula for p s (T)/p, Eq.(36). We assume 
fi — X changes with temperature according to the formula 

2 



H{T) - £(17) = - f*0U ) k B T x 



rp \ 3/2 1 " 



T. 



1 



(46) 



{g a (x) = ^2 n x n /n a ). Here we assume that the particle interaction U and the particle 
density p of liquid helium 4 are renormalized to T\ = 2.17K, which is an approximation 
that dates back to London. Equation. (36) with Eq.(46) predicts a temperature dependence 
of p s (T)/p. Although it bears a qualitative resemblance to p s (T)/p in Fig. 5, it differs from 
Fig. 5 in that p s (T)/p in Eq.(36) remains very small at 3.7K > T > 2.3K (= T\ + 0.2K), but 
it abruptly increases at 2.3K, and reaches a macroscopic number at T\, thus resembling the 
shape of the letter L. Experimentally, as temperature decreases from 3.7K to T\, p s (T)/p 
gradually increases as in Fig.5. This means that Eq.(36) and the approximation behind it is 
too simple to be compared quantitatively to the real system. With decreasing temperature, 
in addition to the particle with p = 0, other particles having small but finite momenta get 
to contribute to the 1/q 2 divergence of as wen - ( m addition to Eq.(30), a new 

Fp(q, 0) including p ^ also satisfies 1 — UFnjO, 0) = in Eq.(29).) Hence, the total p s {T) 



is a sum of each p s (T) over different momenta 38j. The participation of p ^ particles 
into p s {T) is a physically natural phenomenon. For the repulsive Bose system, particles are 
likely to spread uniformly in coordinate space due to the repulsive force. This feature makes 
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the particles with p ^ behave similarly with other particles, especially with the particle 
having zero momentum. If they move at different velocities along the flow direction, the 
particle density becomes locally high, thus raising the interaction energy. This is the reason 
why many particles with p ^ participate in the superfluid flow even if it is just above T\ . 

(3) Equations (35) and (46) with T on = 3.7 K gives us a rough estimate of U as U ~ 
3.4 x 10 _16 erg. This value is approximately close to U s = h 2 / (ma 2 ) = 3.4 x 10~ 17 erg derived 
from the scattering length a = 0.7 nm, which is based on the sound velocity v s = 220 m/s 
and the Bogoliubov formula v s = (h/m)y4^ma. But this U is somewhat larger than U s . 



VI. DISCUSSION 



A. Superfluidity in the non-dissipative and the dissipative flows 

On the onset mechanism of superfluidity, there are physical differences between the non- 
dissipative and the dissipative flows. 

As an example of the non-dissipative flow, we considered the flow in a rotating bucket in 
Sec. 2, in which the quantity directly indicating the onset of superfluidity is the moment of 
inertia 

u _ jf (i.my (47) 

where If is its classical value. 

(1) I z of a superfluid consists only of a linear term of If, and p s (T) appears as a correction 
to its coefficient. 

(2) The change of \ L ~ X T directly affects I z without bein g en hanced, and therefore just 



above T\, the effect of a small but finite p s (T) on I z is small [10]. 

For the dissipative flow, we considered the capillary flow, in which the quantity directly 
indicating the onset of superfluidity is the shear viscosity. Comparing Eq.(25) to Eq.(47), 
we note the following features in v(T). 

(1) v(T) of a superfluid has a form of an infinite power series of v n , and the influence of 
Bose statistics appears in all coefficients of higher-order terms except for the first-order one. 
This feature does not depend on a particular model of a liquid, but on the general argument. 
(On the other hand, the microscopic derivation of the value of v n depends on the model of 
a liquid, which is a subject of the liquid theory and beyond the scope of this paper.) 
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(2) Because of the emergence of the sharp peak due to Eq.(22) in the dispersion integral, 
a small change of p s (T) is strongly enhanced to an observable change of u(T) in Eq.(25). 

(3) The existence of l/d 2 in front of p s (T)/p in the denominator of Eq.(25) indicates that 
a narrower capillary shows a clearer evidence of a frictionless flow. Similarly, the existence 
of l/u e indicates that the choice of experimental procedure such as the method of applying 
the pressure to both ends of the capillary affects the temperature dependence of u(T) of the 
capillary flow. This means that the dissipative phenomena depend on more variables than 
the non- dissipative ones. 



B. Comparison to thermal conductivity 

There is another type of significant change of conductivity in liquid helium 4, an anoma- 
lous thermal conductivity at T < T\. Under a given temperature gradient VT, a heat flow 
q satisfies q = — kVT, where k is the coefficient of thermal conductivity. In the critical 
region above T\ (T/T\ — 1 < 10~ 3 ), the rapid rise of k is observed, and finally at T = T\, 



k jumps to an at least 10 7 times higher value 39J. (This rise is a subject of the fluctuation 
theory of liquid helium 4 {(J.) Comparing to rj(T) which shows a symptom of its fall in 
T\ < T < 3.7 K, k(T) does not show a symptom of its rise in the same temperature region. 

The coefficient of thermal conductivity has a formally similar structure of correlation 
function to that of the shear viscosity, but they are qualitatively different phenomena. 
While shear viscosity is associated with the transport of momentum (a vector), thermal 
conductivity is associated with that of energy (a scalar). For the vector field (velocity 
field), the direction of vectors has a huge number of possibilities in its spatial distribution. 
Hence, among various flow-velocity fields, we can regard the flow in a rotating bucket as 
a non-dissipative counterpart to the capillary flow. On the other hand, for the scalar field 
(temperature field), the variety of possible spatial distributions is far limited. A heat flow is 
always a dissipative phenomenon, and there is no non-dissipative counterpart. Hence, to the 
onset of the anomalous thermal conductivity, the mechanism that amplifies the small p s {T) 
to an observable change cannot be applied. This formal difference between shear viscosity 
and thermal conductivity is consistent with the experimental difference between rj(T) and 
k{T) at T A < T < 3.7 K. 
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C. Comparison to Fermi liquids 



The fall of the shear viscosity in liquid helium 3 at T c is a parallel phenomenon to that 
in liquid helium 4. The formalism in Sec. 2 and 3 is applicable to liquid helium 3 as well. 
For the behavior above T c , however, there is a striking difference between liquid helium 3 
and 4. The phenomenon occurring in fermions in the vicinity of T c is not a gradual growth 
of the coherent wave function, but a formation of the Cooper pairs by two fermions. (This 
difference evidently appears in the temperature dependence of the specific heat: C(T) of 
liquid helium 3 shows only a sharp peak at T c , and does not show a symptom of its rise above 
T c .) Once the Cooper pairs are formed, they are composite bosons with high density at low 
temperature, and immediately jumps to the superfluid state. Hence, the shear viscosity of 
liquid helium 3 shows an abrupt drop at T c without a gradual fall above T c . 

In electron superconductivity, there is an energy gap due to the formation of Cooper 



pairs. 
w = 



3ence, in its conductivity spectrum a{oj) at T < T c , there is a frequency gap uj g near 



40| . In superfluid helium 4, due to the acoustic phonon, there is no energy gap, which 



is consistent with that o n (uj) in Eq.(43) is a weakly tu-dependent function. 
APPENDIX A: h(t) = h(0)exp(-at) 

The mass of a liquid passing through a capillary per unit time is Q = 2np J Q d rv z (r)dr, 
which is determined by Eq.(3) as 

vrd 4 AP /A , 

q = 1»-l- (A1) 

The mass of a liquid in the reservoir with a radius R, pirR 2 h(t), flows out through a capillary 
at the rate of Q. The pressure at the upper and lower end of the capillary set up vertically 
is 1 + pgh and 1 atm, respectively, thus leading to AP(t) = pgh(t) in Eq.(Al). The level 
h(t) of the liquid in the reservoir decreases according to 

pnR^ = -^^h, (A2) 
dt 8i] L K ' 

thereby leading to h(t) = h(0) exp(— at) with a = pgd A j '(8r]LR 2 ). For liquid helium 3, pjr\ 
in a is 3.6 x 10 3 at T = 1.1K, and 3.8 x 10 3 at T = 3.0K. For liquid helium 4, p/r/ is 
4.2 x 10 3 at T = 3.0K. 
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APPENDIX B: CONDUCTIVITY SPECTRUM a(co) 



The Stokes equation under the oscillatory pressure gradient APe wt jL is written in the 
cylindrical polar coordinate as follows 

dv 



dt 



The velocity has the following form 



d d \ APe iu}t 
dr 2 rdr J pL 



(Bl) 



APe iuJt 

under the boundary condition v(d,t) = 0. In Eq.(Bl), Av(r,t) satisfies 



(B2) 



dAv(r, t) 
dt 



= MM), 



(B3) 



and therefore Av(r, t) has a solution written in terms of the Bessel function J (i\r) with 
A = (1 + i)y/u/{2v). Hence, 



v(r,t) 



APe 



iu)t 



icopL 



Jp(iAr) 
Jo(i\d) J 



(B4) 



Using Eq.(B4) at r = 0, the conductivity spectrum a{u), which is defined by Eq.(6) as 
pv(0,t) = a(co)d 2 APe lujt /L, is given by 



a(co) 



icod 2 



1 - 



CO 



J (id[l+i] A /_^ 




(B5) 



The real part of Eq.(B5) 



Rea(co) = 



^2 



n=0 



n\) 2 \ Av 



cod' 



2\ n 



COS • 



cud 2 



oo 



n=0 



;2 \ ™ 

a; a \ n7r 
n!) 2 V W C ° S ^" 



OO j 



n=0 



n!) 2 \ Av 



cod' 



2\ n 



sin 



2 ' 



(B6) 



gives a curve of cr(c<j) in Fig. 3(a). (Re cr(0) in Eq.(B5) agrees with p/(Arj).) Furthermore, 
it determines the conserved quantity f(d) oc d~ 2 in Eq.(7). (f in Eq.(B5) disappears in 
/ a{uo)duo.) 
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APPENDIX C: FLOW IN A ROTATING BUCKET 



The hamiltonian of a liquid in a coordinate system rotating with the container is H—fl-L, 
where L is the total angular momentum. The perturbation H ex = — fl ■ L is cast in the 
form — ^j(f2 x r) ■ p , in which x r = «d(r) serves as an external field. The rotation is 
equivalent to the application of the external field. We define the mass-flow density j{r) in 
the rotation, and express the perturbation H ex as 

-Sl-L = -J v d (r)-j(r)d 3 x. (CI) 

Because of divv d (r) = 0, v d (r) in Eq.(Cl) acts as a transverse- vector probe to the excitation 
of bosons. This fact allows us a formal analogy that the response of the system to v d (r) is 
analogous to the response of the charged Bose system to the vector potential A(r) in the 
Coulomb gauge. Hence, in momentum space, j(r) in Eq.(Cl) has the following microscopic 
form similar to that in the charged Bose system 

U^r) = J2{p + l) $ !V^ r . (C2) 

p,n ^ 

(H—l and r = it). 



APPENDIX D: MAXWELL'S RELATION 



Consider the shear deformation of a solid and of a liquid. In a solid, shear stress F xy is 
proportional to a shear angle as F xy = G<f>, where G is the modulus of rigidity. The value 
of G is determined by dynamical processes in which vacancies in a solid move to neighboring 
positions over the energy barriers. As increases, F xy increases as follows, 

d J3L = G d ±. (Dl) 
dt dt v 1 



In a liquid, the situation is different as follows. Figure 6 represents a small portion of the 
flow in Fig.2, in which two particles 1 and 2, each of which starts at (x,y) and (x,y + Ay) 
simultaneously, moves along the x-direction. In a liquid, the relative position of particles is 
not rigid, but with the flow motion it changes so as to reduce the shear stress F xy between 
adjacent layers. Presumably, the larger shear stress induces the faster structural relaxation, 
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FIG. 6: In a liquid flowing along the x-direction, owing to the velocity gradient along the y- 
direction, a small rectangular part of a liquid is deformed to a parallelogram. 

that is, the rate of such relaxations is proportional to the magnitude of F xy . Hence, instead 
of Eq.(Dl), one obtains 

'ffiry _ _ ^ x v (D2) 

dt dt Tgt' 

where r st is a relaxation time. In the stationary flow after relaxation, F xy remains constant, 
and one obtains 

G d 4 = ^- (D3) 
dt r st 

Assume that there is velocity gradient v x (y) along y direction in Fig.6. After At has passed, 
they (1' and 2') are at a distance of Av x At along the x-direction. The shear angle increases 
from zero to A<f>, which satisfies Av x At = AyA(ft as depicted in Fig.6. Hence, we obtain 

^ = d ± (D4) 
dy dt' K J 

Substituting Eq.(D4) into Eq.(D3), and comparing it with Eq.(l), one obtains r\ = Gr st 
(Maxwell's relation). 
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